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Abstract 

The problem of parameterizing the interactions of disparate scales in fluid flows is ad- 
dressed by considering a property of two-dimensional incompressible turbulence. The prop- 
erty we consider is selective decay, in which a Casimir of the ideal formulation (enstrophy in 
2D flows) decays in time, while the energy stays essentially constant. This paper introduces a 
mechanism that produces selective decay by enforcing Casimir dissipation in fluid dynamics. 
This mechanism turns out to be related in certain cases to the numerical method of antici- 
pated vorticity discussed in Sadourny and Basdevant [1981, 1985]. Several examples are given 
and a general theory of selective decay is developed that uses the Lie-Poisson structure of the 
ideal theory. A scale-selection operator allows the resulting modifications of the fluid motion 
equations to be interpreted in several examples as parameterizing the nonlinear, dynamical 
interactions between disparate scales. The type of modified fluid equation systems derived 
here may be useful in turbulent geophysical flows where it is computationally prohibitive to 
rely on the slower, indirect effects of a realistic viscosity, such as in large-scale, coherent, 
oceanic flows interacting with much smaller eddies. 
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1 Introduction 

1.1 Symmetries and selective decay 

From the viewpoint of Noether's theorem, energy is conserved in ideal fluid dynamics because of 
the time-translation symmetry of the Lagrangian in Hamilton's principle for ideal fluid motion. A 
second type of fluid conservation law arises via Noether's theorem because of relabelling symmetry 
of the Lagrangian. Relabelling symmetries smoothly transform the labels of the fluid parcels 
without changing the Eulerian quantities on which Hamilton's principle for ideal fluids depends. 
The conservation laws associated with relabelling symmetries are called Casimirs, because in the 
Hamiltonian formulation of ideal fluid dynamics in the Eulerian representation their Lie-Poisson 
brackets with any other functionals vanish identically. Thus, the Casimirs arise from a geometric 
symmetry of the Eulerian representation of ideal fluid dynamics. This relabelling symmetry is also 
responsible for Kelvin's circulation theorem in ideal fluid dynamics, which immediately leads to the 
conservation of the Casimirs for 2D ideal incompressible flow, see Holm, Marsden and Ratiu [1998]; 
Cotter and Holm [2012]. The Kelvin circulation theorem and the associated Casimir conservation 
laws are kinematic, because they hold for any choice of Hamiltonian in the Eulerian representation. 
Energy conservation is dynamic. 

The two types of conservation laws for ideal fluids (energy and Casimirs) have been used 
together previously as the basis for the energy- Casimir method of nonlinear stability analysis of 
equilibrium flows HMRW [1985]. This method is based on the concept of dynamical balance in 
ideal fluids. In particular, the equilibrium solutions of ideal fluid dynamics are characterized by 
such a balance. Namely, the critical points 8{h + AC) = of any linear combination h + AC of 
the energy h, a Casimir C and constant A occur for equilibrium solutions of Euler's equations 
for ideal fluids. An initial disturbance of this balance may be either stable or unstable in the 
Lyapunov sense. Sufficient conditions for linear Lyapunov stability of an equilibrium flow may 
be determined by taking the second variation of the energy- Casimir sum, 5 2 (h + AC)| e , where 
subscript e means evaluated at the equilibrium flow. This second variation is conserved by the 
linearized equations. This means that linear Lyapunov stability is implied for those flow equilibria 
for which the second variation S 2 (h + AC)| e comprises a norm in the space of perturbations. Thus, 
in the ideal Hamiltonian description of fluids, the interplay of energy and Casimirs determines 
equilibrium solutions and governs their Lyapunov stability. For more details, and discussions of 
the extension of the energy- Casimir method to nonlinear Lyapunov stability analysis, see HMRW 
[1985] and references therein, particularly Arnold [1965a,b] for the case of 2D ideal incompressible 
fluid flows. 

The selective decay hypothesis in 2D turbulence. Already in Fjortoft [1953] it was shown 
that, in contrast to 3D incompressible flow, the two energy- Casimir constants of motion imply that 
any transfer of kinetic energy to higher wavenumbers (k) in the 2D incompressible flow of a viscous 
fluid must be accompanied by a larger transfer of energy to lower wavenumbers. An explanation for 
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this effect based on statistical mechanics and equilibrium thermodynamics was offered in Kraichnan 
[1967, 1971], which predicted a forward cascade (k~ 3 ) of enstrophy and an inverse cascade (k~ 5 ^ 3 ) 
of kinetic energy in 2D incompressible turbulence, as reviewed in Kraichnan and Montgomery 
[1980]. Investigating these predictions numerically, Matthaeus and Montgomery [1980] performed 
simulations of 2D incompressible turbulence using the Navier-Stokes equations and on this basis 
postulated the selective decay hypothesis. Their selective decay hypothesis was that in the 2D 
incompressible flow of a viscous fluid the squared L 2 norm of vorticity (or enstrophy, the quadratic 
Casimir) should decay in time, while the squared L 2 norm of velocity (kinetic energy) which has 
contributions from both directions in wavenumber should stay essentially constant. The balance 
between the (/c~ 5 / 3 ) and (k~ 3 ) cascades occurs at the length scale at which energy flux passes 
through zero, so that energy may be passed in both directions at the same rate simultaneously. 

Because the balance induced by the energy- Casimir dynamics is so important in the motion 
of fluids, one could imagine designing modified fluid equations that would govern the balance by 
imposing the selective decay of one or the other of these two types of conservation laws. Moreover, 
because this balance in the ideal case stems from the Lie-Poisson bracket in the Hamiltonian 
formulation of ideal fluid dynamics, one could imagine that an effective approach in deriving such 
modified equations would involve using elements of the Lie-Poisson bracket that would select which 
type of conservation law would dominate and thus govern the properties of the solution. 

The two types of ideal fluid constants of motion, energy and Casimirs, typically have quite 
different dependences on spatial gradients of the solutions. Consequently, the interplay between 
them can be interpreted as an interaction between larger and smaller scales (or coherence lengths, 
or spectral wavenumbers). Thus, imposition of a selective decay mechanism for one or the other 
type could in principle control the dominance of one set of scales, or sizes of excitations, over 
another. Hence, modifying the fluid equations to impose selective decay of either the energy or 
a particular Casimir is interesting, because such a modification could control the direction of the 
energy cascade, i.e., whether it is forward (toward smaller, finer scales), or backward (toward 
larger, coarser scales) This means that modifications of the fluid equations that promote selective 
decay could control the dominant sizes that emerge in the flow and support their stability, by 
causing smaller perturbations to decay away. Such modifications would control the rate at which 
the cascade occurs, in either direction, and thus would control the type of solution that emerges 
from a given initial condition. 

Various types of modifications of the Poisson bracket for Hamiltonian systems have previously 
been proposed that produce energy dissipation. These proposed modifications include: (i) adding a 
symmetric bilinear form, Kaufman [1984]; Morrison [1984]; Grmela [1984]; Ottinger [2005]), and (ii) 
application of a double bracket, Brockett [1991]; Bloch et al. [1996]; Holm, Putkaradze and Tronci 
[2008]; Brody, Ellis and Holm [2008]. See also Vallis, Carnevale, and Young [1989], in which a 
modification of the transport velocity was used to impose energy dissipation with fixed Casimirs. 

Here we are not considering energy dissipation. Instead, we are considering a type of selective 
decay of ideal Casimir constants of motion at fixed energy, as was discussed for 2D turbulence 
cascades in Matthaeus and Montgomery [1980]. In the present case, the selective decay is en- 
forced by a modification of the vorticity equation that is based on the Lie-Poisson structure of its 
Hamiltonian formulation in the ideal case Arnold [1966, 1969, 1978]; Holm, Marsden and Ratiu 
[1998]. 

Thus, the aim of this paper is to impose selective decay by using the Lie-Poisson structure of the 
ideal theory, and interpret the resulting modifications of the equations as a means of dynamically 
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and nonlinearly parameterizing the interactions between disparate scales directly. This type of 
modification may be useful in situations where it is computationally prohibitive to rely on the 
slower, indirect effects of viscosity. Moreover, these ideas should apply not only to 2D flow, but 
also to any 3D flow that is sufficiently anisotropic, particularly in geophysical fluid flows in thin 
domains. 



Plan of the paper. The remainder of this introductory section provides more background about 
the two types of conservation laws in fluids and defines our selective decay mechanism based on 
Casimir dissipation in the example of 2D incompressible flows. Remarkably, in that 2D case, the 
mechanism for selective decay by Casimir dissipation in 2D incompressible fluid dynamics turns 
out to be related to the numerical method of anticipated vorticity in Sadourny and Basdevant 
[1981, 1985]. Section 2 discusses the general theory of selective decay by Casimir dissipation in the 
Lie algebraic context that underlies the Lie-Poisson Hamiltonian formulation of ideal fluid dynam- 
ics, as explained in, e.g., Holm, Marsden and Ratiu [1998]. In particular, we develop the Kelvin 
circulation theorem and Lagrange-d'Alembert variational principle for Casimir dissipation. In the 
Lagrange-d'Alembert formulation, the modification of the motion equation to impose selective de- 
cay is seen as an energy- conserving constraint force. Section 3 extends the Casimir dissipation 
theory to include fluids that possess advected quantities such as heat, mass, buoyancy, magnetic 
field, etc., by using the standard method of Lie-Poisson brackets for semidirect-product actions of 
Lie groups on vector spaces reviewed in HMRW [1985]. Here, the main examples are the rotating 
shallow water equations and the 3D Boussinesq equations for rotating stratified incompressible 
fluid flows. Finally, section 4 summarizes our conclusions and suggests some possible next direc- 
tions for applications of the selective decay modifications of fluid equations, treated as dynamical 
parameterizations of the interactions between disparate scales. 



1.2 2D flows 



An interesting feature of ideal incompressible 2D fluid flows is that they have two types of con- 
servation laws which arise from their Hamiltonian formulation in terms of a Lie-Poisson bracket 
{•,•}, Arnold [1966, 1969, 1978], as 

df(u 



dt 



{f,h}{uj) = luj 



~5f 8h 






5h~ 


5u ' 5u 


Sou' 


5u 



where u is the vorticity of the flow, the square bracket 
forms as 



dxdy , (1.1) 
is the 2D Jacobian, written in various 



[f, h] = J(f, h) = f x h y - h x f y = z • V/ x Vh = z • curl (/z) x curl (hz) 



(1.2) 



and the angle bracket ( • , • ) in (1.1) is the L 2 pairing in the domain T> of the (x,y) plane. For 
convenience, we shall take the domain T> to be periodic, so we need not worry about boundary 
terms arising from integrations by parts. 



Energy. The first type of conservation law for such 2D ideal fluid flows arises from the antisym- 
metry of their Lie-Poisson bracket. Namely, antisymmetry of the bracket implies conservation of 
the Hamiltonian function h 

^ = {M>M = o, 
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for any given choice of h. 

Casimirs. The second type of conservation law arises because their Lie-Poisson bracket has a 
kernel (i.e., is degenerate), which means there exist functions C(u) for which 

^ = {^}H = 0, (1.3) 

for any Hamiltonian h(u). Functions that satisfy this relation for any Hamiltonian are called 
Casimir functions. (Lie called them distinguished functions, according to Olver | 0].) 

Theorem 1.1 (Arnold [1966, 1969, 1978]). The Casimirs for the Lie-Poisson bracket (1.1) in the 
Hamiltonian formulation of 2D incompressible ideal fluid motion are 



C$ = / $(u) dxdy, 
Jv 

for any smooth function $. 

Proof. The well-known proof of this statement follows from an identity for the (x, y) Jacobian 
bracket [ • , • ] that arises from integration by parts as 



a[b, c] dxdy = / b[c, a] dxdy = / c[a, b] dxdy . (1.4) 
v Jv Jv 

This identity holds for any smooth functions a,b,c G J-(T>) that vanish on the boundary of the 
domain T>, and also when the domain is periodic. Thus, 

{C$, h}(uj) = / oj $'(cA — dxdy = / — [u, &(u})] dxdy = , 
Jv L " u \ Jv duj 

for any Hamiltonian h and function □ 

Enstrophy. Among the Casimirs C$(a;) for 2D ideal fluids is the famous enstrophy, given by 

1 If 

C 2 ■= ^IMl^v) = 2 / u 2 dxdy, (1.5) 



v 



1, ,2 



which is the Casimir of the Lie-Poisson bracket (1.1) for the case that $(w) = 
Remark 1.2. Any smooth function of a Casimir (such as its square) is also a Casimir. 

Turbulence and selective decay. When viscosity is added, the Navier-Stokes equations for 
2D incompressible turbulence result from the equations above. These equations are dissipative 
and, in the absence of forcing, both types of conserved quantities in the ideal case would decay 
in time and eventually vanish. Remarkably, in 2D incompressible turbulence these two types of 
conservation laws are found numerically to decay at different rates and this difference is found to 
have an important emergent effect on the spectral properties and statistics of 2D incompressible 
turbulence, Kraichnan [1967, 1971]. In particular, this feature of selective decay in 2D turbulence 
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leads to an inverse cascade of energy toward larger scales and a forward cascade of enstrophy toward 
smaller scales where its dissipation occurs, Matthaeus and Montgomery [1980]. This emergent 
effect of viscosity in 2D incompressible turbulence has had a long history of investigation and is 
thought to be important in many applications, including climate modelling. Recent work tends to 
understand this effect as a mechanism for parameterizing the interactions between disparate scales, 
for example, between large coherent oceanic flows and much smaller eddies, Marshall and Adcroft 
[2010]. For other recent discussions of these ideas in modelling disparate scale interactions in 
geophysical and astrophysical turbulence, see Mininni, Pouquet and Sullivan [2008]. 

Question. Given that the two types of conservation laws in ideal 2D flow are both properties of 
the Lie-Poisson bracket for such flows, and their dissipation has a profound effect of the properties 
of the flow, it is natural to ask whether one may use the Lie-Poisson bracket of ideal 2D flow to 
impose a process of selective decay. 



1.3 Casimir dissipation 



The Lie-Poisson bracket may be used to introduce a type of "nonlinear viscosity" that preserves 
the energy, but dissipates a given Casimir. This may be accomplished naturally by modifying the 
Hamiltonian formulation to introduce a quadratic Lie-Poisson bracket structure, as follows, 



df(u) 
dt 



h 

IV 


~s£ 


5h~ 


dxdy — 9 

Jv 


'S£ 


5h~ 


L 


~8C 


5h~ 


5u' 


5u 


Sco' 


5uj 


5u ' 


5u 



dxdy . 



;i.6) 



where 9 is a given constant and L is an arbitrary positive self-adjoint linear differential operator. 
The differential operator L allows a degree of scale dependence in the Casimir dissipation approach. 
For example, choosing L = (1 — a 2 A) s would define a Sobolev H S (T>) inner product with length 
scale a. One immediately sees from equation (1.6) that 



dh(ui) 
dt 



and 



dC{uo) 


—J 


~5C 


5h~ 


L 


~5C 


5h~ 


dt 


Jv 


5uj ' 


5u 




5u ' 


5u 



dxdy =: — 9 



5C_ 6h 

5u ' 5<jj 



The latter equation recovers the Hamiltonian case when 9 = 0. If a given Casimir C(u>) is not 
sign-definite, one may take its square, which is still a Casimir, and find that 



l dC(u) 2 
2 di 



9C{uj) 





~5C 5h~ 




~5C 8h 


t 


Sco ' 5u 


L 





dxdy = -6C{u;f 



5C 5h 
5u ' 5uj 



;i-7) 



for co verifying the equations (1.6) associated to the Casimir function |C(co>) 2 . These calculations 
have demonstrated the following. 



Proposition 1.3. For 9 > in equation (1.6), the energy Hamiltonian h(co) is conserved and 
squared Casimirs C(u) 2 decay exponentially at a rate proportional to the square of the (possibly 
degenerate) norm induced by the L-pairing with linear operator L. 

Remark 1.4. L is called a scale-selective operator, because choosing it to contain higher spatial 
derivatives emphasizes the higher wavenumbers in the L-norm in equation (1.7). This, in turn, 
determines the Casimir dissipation rate. Thus, the choice of the operator L establishes the range 
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of sizes of spatial scales which most contribute to the modifications of the equations due to Casimir 
dissipation. In particular, if one were to take 

L = (1 - a 2 A) s , so that 

then scale sizes greater or less than the length a would tend to behave differently. 





~5C 


Sh' 




z 




~SC Sh' 






Sco ' 


Sco 




L 




Sco ' Sco 





Double Lie-Poisson bracket vorticity dynamics (1.6) in 2D. The 2D vorticity dynamics 
generated by the quadratic Lie-Poisson bracket structure in equation (1.6) is found by choosing 
f(co) = co above, for which 



dco 

~dt = 


'Sh 

Sco 


co 


+ e 


Sh ' 

Sco ' 




'Sh 

Sco 


co 


+ 9L 


'Sh 5C 
Sco' Scu 



Sh SC 
Sco' Sco 



after using the identity in (1.4) obtained from integration by parts. 

Remark 1.5 (Double bracket dissipation structure and "anticipated vorticity"). 



'U 



• When the operator L = Id, equation (1.8) defines a double bracket dissipation structure 
whose formula is reminiscent of the double bracket dissipation introduced in Bloch et al. 
[1996] and Holm, Putkaradze and Tronci [2008] which would yield here the dissipation term 
9 [oj, [a;, J^j]]- Both the similarities and differences between these two approaches can be 
clearly seen from their respective abstract Lie algebraic formulations, see Remark 2.6 later. 

• When one sets Sh/Suo = ip, with stream function ip(x,y), and chooses the Casimir to be the 
enstrophy C — 5 j v 0J 2 dxdy, equation (1.8) recovers the anticipated vorticity model (AVM) of 
Sadourny and Basdevant [1981, 1985] for an appropriate choice of L. See also Vallis and Hua 
[1988]. Namely, in terms of the stream function representation of the 2D vector velocity 
u = curl (ipz) = — z x with vorticity 00 z = curl curl (ipz) = —Aifj z, equation (1.8) above 
with L = Id and Sh/Su = ip becomes 

^ = [V , co] + 9 [V , L [V, co]] = {ij , co + 9L ty, to]] 

= - u ■ Vw + 9 u • VL(u • Vw) = - u • V(w - 6 L(u ■ Vco)) . 

The word "anticipated" was used in Sadourny and Basdevant [1981, 1985] when naming 
AVM, because if one were to take 9 ~ At, where At is the time step of the numerical method, 
the term proportional to 9 in the last parenthesis of (1.9) would approximate the vorticity at 
t + At for the pure Hamiltonian vorticity evolution at linear order in At. The AVM approach 
has been very intriguing in the modelling of 2D incompressible turbulence and much has been 
written about it in that literature. These ideas and particularly the choice of the time-scale 
parameter 9 have also been discussed recently in the context of shallow water modelling, see, 
e.g., Chen, Gunzburger, and Ringler [2011a,b]. The other scale- selective AVM models con- 
sidered in Sadourny and Basdevant [1981, 1985] can be obtained by making an appropriate 
choice for the differential operator L in equation (1.6), or (1.8). Sadourny and Basdevant 
[1981, 1985] take L = A 8 , corresponding to the Sobolev space H 8 . 
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Remark 1.6 (Quasigeostrophy). The flexibility of the Casimir dissipation approach can be shown 
by passing to the case of 2D rotating quasigeostrophic fluids (QG), for which the vorticity uj is 
replaced by potential vorticity q satisfying the relation 5h/5q = xjj and equation (1.9) becomes 

^ = -u- V(q-6L(u- Vg)), u = -zxVV>, q = -Ai/j + Tip + f , 
ot 

where the constants T and / denote the square of the inverse Rossby radius and the rotation 
frequency, respectively. 



Double Lie-Poisson bracket vorticity dynamics in 3D. Our goal now is to write the 
quadratic Lie-Poisson bracket structure (1.6) in any number of dimensions and illustrate its effects 
in several hopefully illuminating examples. We begin by using (1.2) in 2D to extend the quadratic 
vorticity bracket expression (1.6) from 2D to 3D, as follows, 



dt 



{/,*}(«) 



J/ ,6h j3 
uj ■ curl - — x curl - — d x 

OUJ OUJ 



— v I curl - — x curl - — • L curl - — x curl - — d" 

J v OUJ OUJ \ OUJ OUJ J 



;i.io) 



X . 



The first (skewsymmetric, Lie-Poisson) summand in this modified vorticity bracket may be written 
in terms of velocity equivalent ly as 



{/, h}xst(u) 



J/ Jh j3 
uj ■ curl - — x curl - — a x 

OUJ OUJ 



u ■ curl 



v 



Sf 5h 

ou Oil 



d 3 x 



- / u 

IV 



5£ 5h 
5u ' 5u 



d A x . 



where now in 3D the square brackets [ • , • ] represent Lie brackets of divergence-free vector fields, 
and we have used the identity [u, v] = — curl(u x v) for V • u = = V • v. 

The second (inner product) summand in the vorticity bracket (1.10) depends on the Casimir 
C and may be written in terms of velocity equivalently as 



r , , / n nf Jf Jh T ( JC ,Sh\. 
\ t , n\ond{xi) — —vl curl - — x curl - — • L curl - — x curl - — d" 
u ' il y ' J v 5uj 5uj \ 5uj 5uj) 



x 



-9 



v 



5f 5h 



/SC 5h 



■f- x — ) • L I — x — I d 6 x 
ou ou 



ou ou 



V Su 5u 

curl _1 Lcurl _1 curl ( — — x — — 



5f 5h 



1-1 T 1-1 



SC Sh 



\ 5u 



Su 



d s x 
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Here, the positive operator A := curl x Lcurl 1 defines a symmetric bilinear form 7 on velocities. 
That is, on the Lie algebra of divergence-free vector fields we have 7 : X x X — > R 



u,v) := / U'AvA:=(u, v b ) 
Jv 



7( 



The second equality here defines a pairing ( • , • ) : X x X* — > R, with the flat-operator b : X — >• X* 
given by = (Av) ■ (fx. 

Having investigated the structure of the modified vorticity bracket for vorticity uj in equation 
(1.10), we may now write the modified equation for the vorticity, as 



£ = curl ( curl^ x w V cur. fcurl^ x L four. £ x cur. «L 
at \ ouj J \ ouj \ ouj ouj 



curl 



5h 



5u> 



curl- — x ( uj — 9 L ( curl-: — x curl - — ) 



,8C 



8h\ 



\ 5u) 



5ll> J J 



Xll) 



Conserved quantities in 3D. In 3D the fluid kinetic energy is expressed in terms of its vorticity 

as 

h(uj) = - [ \u\ 2 d 3 x = - [ u> ■ (-A _1 o;) d 3 x (1.12) 
and the only Casimir in the Lie-Poisson Hamiltonian formulation is the helicity 

C(u) = -/ uj-ud 3 x = - uj- curl" 1 uj d 3 x , (1.13) 
2 Jp 2 J v 

in which we may take curl -1 u; = u to be divergence- free. The helicity represents the total linkage 
of the lines of vorticity with itself. It is preserved by the Euler fluid equations, but is dissipated 
by the viscous Navier-Stokes equations. For more details and discussions of the nature of fluid 
helicity, see Arnold and Khesin [1998]. For the choices of h(uj) and C(uj) above, we have 

5h _5h . 5C _8C . 

— = curl- — = u and — — = curl- — = curlu = uj . (1-14) 

OU OUJ OU OUJ 

Consequently, the vorticity dynamics equation (1.11) for these choices becomes 



duj 

— - = curl (u x (uj — 9 L (uj x u))) 

dt 

u , uj — 9 curl (A [u , uj] ) 



1.15) 



and when the operator L = Id, we have A = —A 1 (the inverse Laplacian) and curl A = curl . 
In velocity form, these equations read 

d t u + (u>- 9L(uj x u)) x u = -Vp (1.16) 

Remark 1.7 (Stationary solutions and generalized Beltrami flows and Lamb surfaces). 
Stationary solutions of (1.16) satisfy 

(uj — 9L(uj x u)) x u = — Vp . 
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This stationary flow condition means that either the vectors of velocity u and modified vorticity 
u> — 9L{uj x u) are parallel, a subcase of which is u> x u = (Beltrami flows), or that the vectors 
u and us — 6L(uj x u) are both tangent to the same pressure surface p(x) = const, viz 

u-Vp = and (u> - 9L(u: x u)) ■ Vp = . (1.17) 

For 9 = 0, a surface defined by p(x) = const for steady flows would be called a Lamb surface, 
see, e.g., Holm [2011]; so surfaces satisfying (1.17) may be called generalized Lamb surfaces. The 
geometry of generalized Beltrami flows and Lamb surfaces may be an interesting research topic. 

Remark 1.8. Equation (1.15) is the 3D version of the anticipated vorticity equation (1.9) in 2D 
with L replaced by curl A = Lcurl -1 . Again, the constant 9 ~ At has units of time. A "scale- 
aware" version of this anticipated vorticity method in 3D may be obtained by selecting the value 
of the time-scale parameter 9 based on local mean properties of the solution, instead of the time 
step of the numerical method. However, However, then the issue arises of whether the time-scale 
parameter 9 should be part of the symmetric operator L and, thus, whether its spatial dependence 
should play a role. This is beyond the scope of the present work. For a discussion of a numerical 
algorithm with scale awareness for 2D multiresolution grids, see Chen, Gunzburger, and Ringler 
[2011b]. 



Anticipated Kelvin circulation. Inverting the curl in equation (1.15) yields 

(d t + C u )(u ■ dx) = —dp + 9 C u ( curl _1 L(u; x u) • dx) 

= -(ip + e£ M (curr 1 Lcurr 1 (curl(w x u)) • dx) (1.18) 
= —dp + 9 C u (A[u, u>] ■ <ix) , 

where, we have substituted A := curl _1 L curl -1 and used [u, u>] = — curl(u x u>). The symbol C u 
denotes the Lie derivative with respect to the velocity vector field u — u • V, expressible when 
applied to 1-forms in 3D as 

C u {y ■ dx.) = — (u x curl v) ■ dx + d(u ■ v) . (1-19) 

For an introduction to the use of Lie derivatives in fluid mechanics, see Holm [2011]. For a more 
advanced discussion, see Holm, Marsden and Ratiu [1998]. 

Kevin's circulation theorem now becomes, cf. equation (1.15) 

— (f u • dx = <L (d t + C u )(u ■ dx) = —9 (j) u x curl (A[u, u>}) ■ dx 
dtj c{u) / c(u ) / c(u ) 

= —9 (p u x L(uj x u) • dx . 

Jc(u) 

The last term is the source of circulation due to Casimir dissipation. 



Variational-derivative expressions. Recall that, in terms of the velocity, we have the varia- 
tional expressions 

— = u and -r— = u> (C = helicity), 
ou ou 
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so the last line of the calculation in (1.18) may be expressed in terms of variational derivatives as 



(d t + C Sh /su) (u • dx) = -dp + 9 Csh/su 



Sh 5C 
<5u' 8u 



;i.2i) 



This type of expression in terms of variational derivatives will also be obtained by using the general 
theory developed in the next section. 

Proposition 1.9 (Squared Casimir dissipation - Helicity). Let the 3D version of equation (1.7) 
for squared Casimir dissipation be applied to the square of the helicity, with 



C(oj) = - I u) ■ curl 1 ujd z x and h{us 



1 



v 



curl 1 uj\ 2 d 3 x = - I IuPgPx. 



V 



12 j3. 



V 



Then the squared helicity Casimir dissipates according to the equation 

\dC(uS) 2 
2 dt " 



-9C(uj) 2 \\[u,u] 



1.221 



|A -0C(u,y\\VLXU,\\' L . 

Proof. The proof is a direct calculation, starting from equation (1.10). 



ldC(w) 



2 dt 



5C 



8C (a?) / curl - — x curl - — • L I curl - — x curl 



v 



5uj 



5C 



5h 



-9C(u>) 2 [ curl" 1 
Jv 

-0CV 2 



5a; 



curl - — , curl - — 

OUJ OUJ 



5u) 5u> 
■ L curl -1 



d s x 



JC Sh 
curl - — , curl - — 

OUJ OUJ 



d A x 



V 



JC Sh 
curl - — , curl - — 

OUJ OUJ 



A 



JC Jh 
curl - — , curl - — 

OUJ OUJ 



d A x 



1.23) 



-9C(ujY [uj,u] ■ A[uj,u] d A x with A := curl ^curl 



i-i 



v 



= -9C(uj) 2 1 ( k [uj,u},[uj,u] 
=: -6C(uj) 2 \\ [«,u] 

A shorter proof may be obtained by replacing the functional derivatives directly in the first line 
to get the equivalent result, 

ldC(uj) 2 2 2 

and then substituting the interesting identity that relates the L- and A-norms, 



\UJ X u 



[^> U J Ha- 



1.24) 
□ 



2 Selective decay by Casimir dissipation: General theory 

It is clear that the modification process in the previous section produces the anticipated equations 
in the cases discussed so far, but despite appearances in those examples, this Casimir dissipation 
process involves considerably more than merely iterating the evolution operator. 
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2.1 The modified Lie-Poisson (LP) framework 

To describe the modified Lie-Poisson structure for the general theory, we fix a Lie algebra g with 
Lie brackets denoted by [ ■ , ■ ], and consider a space g* in weak nondegenerate duality with g, that 
is, there exists a pairing (•,•): g* x g — >• R, such that for any £ G g, the condition (/i, £) = 0, for 
all n G g* implies £ = and, similarly, for any /i G 0*, the condition (/x, £) = for all £ G implies 
/i = 0. Recall that g* carries a natural Poisson structure, called the Lie-Poisson structure, given 
by 

~Sf Sh" 

(Marsden and Ratiu [1994]) where /, g G J^id*) are real valued functions defined on g* and Sf/S/j, G 
g denotes the functional derivative, defined through the duality pairing (•,•), by 



{M}+0) = (n 



(2.1) 



8f 6 \± 
5ji ' / cfe 



f(n + e5fx). 



£=0 



This Poisson bracket is obtained by reduction of the canonical Poisson structure on the phase 
space T*G of the Lie group G with Lie algebra g. The symmetry underlying this reduction is given 
by right translation by G on T*G. In the case of ideal fluid motion, this symmetry corresponds to 
relabelling symmetry of the Lagrangian in Hamilton's principle. 

Recall that a function C : g* — > R is a Casimir function for the Lie-Poisson structure (2.1) if it 
verifies {C, /}+ = for all functions / G ^(g*) or, equivalently ad*sc \x = 0, for all \i G g*, where 

5/J, 

ad^ : g* — > g* is the coadjoint operator defined by (ad^ yu,, 77) = (/i, [£,^]). 

Below, we will denote by 7 M be a (possibly /^-dependent) symmetric bilinear form 7^ : gxg — > R. 
This form is said to be positive if 

7m(£,£)>0, foralUGg. 

Definition 2.1. Given a Casimir function C(n), a positive symmetric bilinear form 7 M , and a real 
number 6 > 0, we consider the following modification of the LP (Lie-Poisson) equation to produce 
the Casimir dissipative LP equation: 



dtfi + adlh /i = 9 &d*i 



6 j.i 



5C 5h 

5jJL ' 5jjL 



(2.2) 



where b : g — > g* is the flat operator associated to 7^, that is, the linear form £ b G g* g'^en 6y 
£\v) = 7 M (^> /° r al1 k,V^S- 



Note that the flat operator need not be either injective or surjective. Note also that in the 
equations (2.2) above, the flat operator is evaluated at fi. It is important to observe that the 
modification term depends on both the given Hamiltonian function h and the chosen Casimir C. 
It is convenient to write (2.2) as 

dtn + &d*sh fi = 0, (2.3) 

Sfi 

r 1 b 

for the modified momentum \x := fi + 8 |£ 
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Remark 2.2 (Left-invariant case). Recall that the Lie-Poisson structure (2.1) is associated to right 
G-invariance on T*G. We have made this choice because ideal fluids are naturally right-invariant 
systems in the Eulerian representation. Other systems, such as rigid bodies, are left G-invariant. 

In this case, one obtains the Lie-Poisson brackets {/, g}^(u) = — (^fx, ^ \ and this leads to 

the following change of sign in the Casimir dissipative LP equation (2.2): 

d t fj, — a,d*sh fx = 9 ad^k 



8C 5h 
8 fx ' 5 ix 



(2.4) 



The rigid body example. For the rigid body case itself, one identifies fx in (2.4 ) with II G R 3 , 
the body angular momentum, and 5h/5fx with Q e R 3 , the body angular velocity, while the 
Casimir is C = ^|II| 2 . Then, choosing for 7 the usual inner product on R 3 , (2.4) becomes 



dU 

~~dt 



+ n x n = -en x (n x n) . 



Hence, the Casimir is dissipated by this modified rigid body equation when 9 > 0, since 

d 1 



dtT 



ini 



-e n ■ n x m x ft) 



-9\n x ni 2 



(2.5) 



(2.6) 



Remark 2.3 (Casimir dissipation VS anticipated setting). Clearly, equation (2.2) is not equivalent 
to a simple iteration of the evolution operator, as might be considered in an "anticipated" setting 
in the form 

d t /x + ad^, n = 9a,d*sh (ad^k fx) ■ (2.7) 

Sfx Sfx \ Sfx / 

However, the "anticipated" modification does not guarantee selective decay of Casimirs, dC/dt < 0, 
in general. The dynamics of equations (2.2) for Casimir dissipation and (2.7) for "anticipated" 
motion are only equivalent when 



ad* 



6 1 1- 



5C 5h 

5 Li ' 5 Li 



&d*sh ( adih Li 

5fi 



(2.i 



which can only hold for quadratic Casimirs C. Thus, "anticipated" motion and selective Casimir 
decay can only coincide for certain quadratic Casimirs. 



Lie-Poisson formulation. In Lie-Poisson form, the Casimir dissipation equation (2.2) reads 

df(Lx) /5f 



dt 



Ofl 



$f ,* \ n 6 f 1* 

— , adih fx ) +9 I — , a.dsh 

Ofl iH / \ Ofl *M 



5C 5h 
8fi ' 8fi 



~5f 5h~ 




~6_l 


8h~ 




'SC 8h' 


8fi 


>- 


5fx 


8fi 




8 fx ' 8 fx 



(2.9) 



( 


~8f 8h 




~8C 8h 




8 fx' 8 fx 


1 


8 fx ' 8 fx 
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for an arbitrary function / : Q* — > R. The energy is preserved, since we have 

dh(fi) 



dt 



{h,h} + -0 7 



~5h 


5h~ 




~8C 8h 




5fi 




5fi ' 8fi 



0. 



However, when 9 > the Casimir function C is dissipated since 
dC(n) 



dt 





~5C 


5h 




~5C 


5h~ 


)- 




"5C <S/f 




( 


5fi 


5fi 




5fi 


5fi 




<5/A ' 5/i 





(2.10) 



where ||£||^ := 7 M (£, is the quadratic form (possibly degenerate) associated to the positive bilinear 
form 7 M . 



Left invariant case. In the left invariant case, we have, in comparison with (2.9). 



dt 



{f,h}_-6j 



( 


'Sf 


5h 




~8C 


5h~ 








5fi 




5fi 


5fi 


) 



Kelvin-Noether theorem. The well-known Kelvin circulation theorems for fluids can be seen as 
reformulations of Noether's theorem and, therefore, they have an abstract Lie algebraic formulation 
(the Kelvin-Noether theorems), see Holm, Marsden and Ratiu [1998]. We now examine how the 
dissipation term introduced above modifies the abstract Kelvin circulation theorem. 

In order to formulate the Kelvin-Noether theorem, one has to choose a manifold C on which the 
group G acts on the left and consider a G-equivariant map /C : C — > Q**, i.e. (/C(gc), Ad*_i v) = 
(/C(c), v) , V g G G. Here gc denotes the action of g G G on c G C and Ad* denotes the coadjoint 
action defined by ^Ad*/i, £) = (/iAd 9 £), where /i G g*, (Eg, and Ad 3 is the adjoint action 
of G on g. Given c G C and /i G g*, we will refer to (IC(c),fi) as the Kelvin-Noether quantity 
(Holm, Marsden and Ratiu [1998]). In application to fluids, C is the space of loops in the fluid 
domain and /C is the circulation around this loop, namely 

(/C(c), u ■ dx) := j) u • <ix. 

The Kelvin-Noether theorem for Casimir dissipative LP equations is formulated as follows. 

Proposition 2.4. Fix cq G C and consider a solution /i(£) of the Casimir dissipative LP equation 
(2.2). Let g(t) G G be the curve determined by the equation = gg^ 1 , g(0) = e. Then the time 
derivative of the Kelvin-Noether quantity (/C(g(t)co), /i(t)) associated to this solution is 



d 
dt 



(JC(g(t)c ),fi(t)) 



9 ( JC(g(t)c ), &d*6h 



5C 5h 
5fi ' <5/i 



Proof. The proof is a direct calculation. We have 

' (fc(g(t)co),v(t)) = — (£(co),Ad* (t) /i(t)) = ^/C(c ),Ad* (t) (<9 t /i + adk //) 



9( K(g(t)co),ad*ih. 



5C_ 5h 
5fi ' <5/i 



where, at the first equality we used the G-equivariance of K and at the second equality we used the 
formula ^ Ad*( t ) /i(t) = Ad*^(d t n(t) + ad*^^-! /i(t)), see, e.g., Marsden and Ratiu [1994]. □ 
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Note that g(t) EG is the motion in Lagrangian coordinates associated to the evolution of the 
momentum fi(t) G g* in Eulerian coordinates. The 9 term is an extra source of circulation with a 
double commutator. This term is absent in the ordinary Lie-Poisson case (i.e. 9 = 0) and therefore 
in this case the Kelvin-Noether quantity (/C(g(i)c ), n(t)) is conserved along solutions. 

Remark 2.5 (Three dimensional ideal flows). The Casimir dissipation approach for 3D flows 
developed in §2.2 follows from the present abstract formulation by choosing the Lie algebra g = 
'X-divi'D) of divergence free vector fields on P, the helicity Casimir C, and the positive symmetric 
bilinear form 7(11, v) = L, u • Av<i 3 x, where A = curl - L curP with L an arbitrary self-adjoint 
positive operator. Using v b = Av • <ix, it is readily checked that with these choices, equations 
(2.2) and (2.9) yield (1.10) and (1.15). Note that Proposition 2.4 applied to this case yields the 
circulation theorem 

■j- /" u • dx = 9 [ C u (A [u, uj] ■ <ix) = — 9 f u x L{uj x u) -dx, 

which consistently recovers our previous expression for Kelvin's circulation theorem with Casimir 
dissipation obtained (1.20) by a direct calculation. 

Remark 2.6 (Casimir dissipation VS double bracket dissipation). The Casimir dissipative LP 
setting considered here is essentially different from the double bracket dissipation setting. Indeed, 
double bracket equations dissipate energy while they preserve the Casimirs. This is exactly the 
opposite with what happens in the present setting. There are however apparent similarities in 
the Lie algebraic formulations. Indeed, on general Lie algebras the double bracket dissipation 
equations can be written as 

^ = {/, h}+(ii) - fadk I*, ad*,, /x) , (2.11) 
at v <*M <sm / 

(compare with equation (2.9)) where 7 is a inner product on g, 7* is the inner product induced 
on g*, and k : g* — » M is a given function. One readily checks that Casimir are preserved while, 
in the special case k = h, the energy dissipates. In that case, the equation of motion arising from 
the double bracket in (2.11) is given by 

d t jj, + adk ix = 9 ad*, .(t/i, (2-12) 

see Bloch et al. [1996]; Holm, Putkaradze and Tronci [2008], where (j : g* — > g is the sharp oper- 
ator associated to 7. Formula (2.12) for double bracket dissipation may be compared with the 
corresponding equation for Casimir dissipation in (2.2) 



d t jJL + &d*sh [i = 9 adlih 



SC 5h 



(2.13) 



(S /.i 

and also with the corresponding equation for "anticipated motion" in (2.7) 

dtii + adk = 9 adtih (adk fi) . (2-14) 
The last three of these formulas all coincide with (1.9) for the case of 2D incompressible flow. 
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2.2 Lagrange-d'Alembert variational principle. 



We now explain how the Casimir dissipative LP equations can be obtained via a variational prin- 
ciple. Consider the Lagrangian £ : q — y M. related to h via the Legendre transform, that is, we 
have 

5£ 

where we assumed that the second relation yields a bijective correspondence between £ and fi. In 
terms of £, equation (2.2) for Casimir dissipation reads 



dtix + adg fi = 9 ad^ 



6C " 



5£ 



(2.15) 



d t — + ad: — = 



When 9 = we recover the Euler-Poincare equations 

5£ 

and it is well-known that these equations can be obtained via the constrained variational principle 
Arnold and Khesin [1998]; Holm, Marsden and Ratiu [1998] 

f T 

8 / £{£)dt = 0, for variations 6£ = d t ( — [£, Q, 
Jo 

where ( G is an arbitrary curve vanishing at t = 0, T. 

External forces /(£) 6 g* can be included in the Euler-Poincare equations by using the reduced 
Lagrange-d'Alembert principle Bloch [2004] 



5 



[ mdt] + [ (/(£), C) dt = 0, for variations 5£ = - [£, C], (2.16) 
Jo J Jo 

where ( G Q is an arbitrary curve vanishing at t = 0,T. Substituting into (2.16) the formula 
for the force appearing in the Casimir dissipation equations (2.15) recovers equations (2.15), now 
rederived from the Lagrange-d'Alembert variational principle 



£(0dt 



+ 9 



7 



5C_ 

Sfj, 



[£, £] ) dt — 0, for variations 5£ = d t ( — [£, (]. 



Thus, in the Lagrange-d'Alembert formulation, the modification of the motion equation to impose 
selective decay is seen as an energy- conserving constraint force. 

For example, for three dimensional ideal flows, equations (1.15) can be obtained from the 
variational principle 



o Jv 



u\ 2 d 3 xdt + 9 



o Jv 



([(J, u] • A [u, v]) d 3 x dt = 0, for variations 5u = <9 4 v + [u, v]. 



The second term on the left hand side can be also written as 



9 



o Jv 



(ui x u) • L(u x v) d 3 x dt . 



This is the time integral of the rate of virtual work done by the force for a virtual motion v 
5(f o if' 1 of the fluid configuration ip :T> —y D. 
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3 Semidirect product examples 

The Hamiltonian structure of fluids that possess advected quantities such as heat, mass, buoyancy, 
magnetic field, etc., can be understood by using Lie-Poisson brackets for semidirect-product actions 
of Lie groups on vector spaces. 

In this setting, besides the Lie group configuration space G, one needs to include a vector space 
V on which G acts linearly. Its dual vector space V* contains the advected quantities. From 
this, one considers the semidirect product G®V with Lie algebra fj(§)V, and the Hamiltonian 
structure is given by the Lie-Poisson bracket (2.1), written on (g©^)* instead of q*. We refer 
to Holm, Marsden and Ratiu [1998] for a detailed treatment. Given a Hamiltonian function h = 
h(fx,a) on ($©V)* one thus obtains the Lie-Poisson equations 



dtifji, a) + ad*, sh sh\fi = 0, 

\ Sfi 3 8a } 

for fx(t) G Q* and a(t) G V*. More explicitly, making use of the expression of the ad*-operator in 
the semidirect product case, these equations read 



d t fx + &d*sh fx 

8f.L 



5h 
8a 



o a 



0, d t a + a 



5h 
8 fx 



0. 



(3.1; 



(a£,v), and a£ G V* denotes the 



where the operator o : V x V* — > q* is defined by (v O a, £) = 
Lie algebra action of £ G g on a G V*. 

We can easily extend the Casimir dissipation approach of §2 to the semidirect product case. 
Fixing a Casimir function C = C(fi,a) for the Lie-Poisson bracket on the semidirect product and 
a (possibly (fx, a)-dependent) positive symmetric bilinear form 7( M , a ) : x g — > M, we extend the 
semidirect-product Lie-Poisson system (3.1) naturally to allow for Casimir dissipation by setting 



d t fx + &d*sh it + — o a = 6 adk 

<5m oa 



5C 5h 
5 fx ' 5 fx 



dta + a 



5h 
5fx 



0. 



(3.2) 



where b : Q — > Q* is the flat operator associated to 7. In Lie-Poisson form, this equation reads 

df(n,a) 



dt 



{f, h} + (fx, a) -6-1 







5h~ 




~8C 


5h~ 


( 


5 fx 


5 fx 




Sfx 


5 fx 



which shows that the energy is conserved while the Casimir dissipates as in (2.10). 

To formulate the Kelvin- Noether theorem, we consider a G-equivariant map /C : C x V* — > g**, 
i.e. \JC(gc, ag^ 1 ), Ad*_i v) = (JC(c,a),v) , V g G G. Given c G C and a solution fx(t), a(t) of (3.2), 
the associated Kelvin-Noether quantity reads (JC(g(t)c , a(t)), fx(t)) and similar computations as 
in Proposition 2.4 yield the Kevin-Noether theorem 



j t (JC(g(t)c ,a(t)),fx(t)) = { /U (/(/)(,, ait)). adlw. 



5C_ 6h 
5 fx ' 5 fx 



5h 
8a 



o a 



(3.3) 



Remark 3.1 (Natural generalization). Note that the Casimir dissipative LP equation (3.2) is 
obtained from the Lie-Poisson equations on the semidirect product by modifying the momentum 
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11 G g* only, while keeping the quantity a EV* unchanged. From the Lie algebraic point of view, 
however, the direct generalization of (2.2) to semidirect product Lie groups is 



d t (n, a) + ad/ih sh\fi = 9 ad/ 

Sh 5h\ 

\ Sfi ' 6a ) \ Sfi ' 5a J 



5C 5C_\ fbh 5h\ 
5u,' 5a J \5u,' 5a J 



(3.4) 



where the flat operator b : g x V — > g* x V* is associated to a positive symmetric bilinear map 
7(«,a) : (p V) x (g x V) 4 I. This results in a modification of both /i and a as 



/! = /i — 9 



5C 5h 

5jjL ' 5jjL 



Q f5C5h _6h6CY 
\ 5a 5fi 5a 5ji J 



Of course, as before, the modified semidirect-product Lie-Poisson system (3.4) introduces dissi- 
pation of the Casimir C while keeping energy conserved. It recovers (3.2) in the case when the 
bilinear form 7 vanishes on V. 

Lagrange-d'Alembert variational principle for semidirect products. As in (2.15), equa- 
tions (3.2) for Casimir dissipation in the semidirect product case can be expressed from the La- 
grangian i : g x V* — > R associated to h via the Legendre transformation 



h((M,a) = (//,£) - £(£,a), n := 



(3.5) 



When 9 = 0, the constrained variational principle associated to the equations reads 
f T 

5 / a)dt = 0, for variations 5£ = d t ( — [£, C], 5a = —a(, 
Jo 

where ( G g is an arbitrary curve vanishing at t = 0,T, see Holm. Marsdcn and Ratiu [1998]. 
As one may verify directly, when 9^0, the dissipative Casimir force can be included in the 
Euler-Poincare equations by considering the Lagrange-d'Alembert variational principle 



6 



£{$,, a)dt 



+ 9 f 1 
Jo 



5C_ 
5[i 



[£,C] U = o, 



(3.6) 



for variations 5£, = d t ( — [£, (], 5a = —a(. 



3.1 Rotating shallow water (RSW) flows 

On the semidirect product of Diff(D) with functions n G -7 r (^ > ), where Diff(X>) is the group of 
diffeomorphisms of a two dimensional domain V, the Lie-Poisson equations (3.1) become 



m 



m 



6h 



dt[ — )+ curl ( — 1 x — + V + — 



5h m 5h 



V V J V V J ^ m \5m 77 5r] 

The Hamiltonian for the rotating shallow water (RSW) equations is 



Sh 

0, d t rj + div I 77- — 
dm 



0. 



(3.7) 



h(m,r]) = j (^-\m-riR\ 2 + ^g(r]- D) 2 ^j dxdy 



(3i 
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where curlR = 2f2 is the Coriolis parameter and D(x,y) is the mean depth of the bottom to- 
pography in the two dimensional domain with spatial coordinates (x,y). Indeed, inserting this 
Hamiltonian into the Lie-Poisson equation (3.7) yields the rotating shallow water equations for the 
fluid velocity u = m/rj — R and the total depth rj, 

d t u + curl(u + R) x u + V Ql u | 2 + d(v ~ DU = 0, d t n + div(^u) = 0. 

The second term may be written as r]~ l curl(u + R) x nu = qz x rju, with the potential vorticity 
(PV) q given by 

q = r\~ x z ■ curl(m/?7) = r/ _1 z ■ curl(u + R) 

As a consequence of the shallow water equations, the potential vorticity q is advected (conserved 
on fluid parcels), that is 

d t q + u • Vg = . 

Of course, the conservation of PV is not limited to the shallow water equations. In fact, it holds 
for the equations derived from this LP structure for any Hamiltonian. This is because the LP 
structure admits PV in a family of Casimir functions that Poisson commute with every functional 
of the variables (m, 77). 

One class of Casimir functions is given by 

C$(m, 17) = / rj $ (q,r)) dxdy, where q = r]~ x z ■ curl (m/77) , 
Jv 

for any smooth function $. We will compute the Casimir dissipative version of RSW associated 
to the Casimir function 

C ( m ,v) = ^J riq 2 dxdy, (3.9) 

and the inner product 

7jj(u, v) = / 77 (u • v)dxdy . (3.10) 
Jv 

The flat operator associated to this inner product is thus v b = r?v • dx. From the abstract Lie 
algebraic formulation (3.2) we obtain 

<9 t u+curl(u+R+6»[u,zx(Vg)/?7])xu+V Ql u | 2 + d(v ~ D) + 9 u • [u, z x [Vq)/rj\j = 0. (3.11) 

This equation is obtained from the Lie-Poisson equations (3.7), where in the second and third 
terms the expressions m/n have been replaced by m/77, where 



m = m 



Sh 5C 
5m ' 5m 



m + 9 [u, z x (Vq)/r]] b = m + 9r] [u, z x (Vq)/rj\ , 



as dictated by the Lie-Poisson Casimir dissipation approach in formula (3.2). We have used the 
formula 5C/5m = curl (qz) /rj = —z x (Vq)/n for the variational derivative of the Casimir (3.9) 
and kept in mind that the Lie algebra bracket appearing in (3.2) is minus the Lie bracket of vector 
fields. The last term of (3.11) is computed as follows 

5h m 5h , _ „ r „ .„ , . 1 . l9 „ . - , 

7 + -=- = u - (u + R + u,z x CVq) rj }) - - \u\ — u ■ K + g(rj — D) 

5m r] 5r] 2 

= i|u| 2 + g(rj- D) + 9u- [u,z x {Vq)/rj\, 
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where we used the relation u = m/n — R between the velocity and momentum in the rotating 
case. 

This model preserves the energy and dissipates the Casimir \ J v r\ q 2 dxdy as 

^,M = -»/,[»(v,),,.f W ,. 

dt 2 J v J v 



Kelvin circulation theorem. A convenient way to derive the Kelvin circulation theorem (3.3) 
is to rewrite the Lie-Poisson equations (3.7) using the duality pairing with one-forms. We get 



a 



V 



. a\ Sh 
+ £ sh — + d— 
«<* \ rj J on 



0, d t r] + div ( 7] 



Sh 
1 1 5a 







(3.12) 



in which a denotes the 1-form a = m • cfx and Csh/Sa denotes the Lie derivative with respect to 
the vector field 8h/8a, as defined for 3D earlier in (1.19). The Casimir dissipative LP equation is 
obtained by replacing a by a = m • <ix in the second term. Consequently, a direct computation 
yields 



d f a 

dt J c{u) rj 



9 d> Csh 

lc(u) Sa \V 



5h 5C 
Sex' Sa 







c(u) 



jCsh 

5a 



Sh 5C 
5a' Sa 



consistently with (3.3), where we recall that the Lie algebraic bracket is minus the Lie bracket of 
vector fields. For RSW, the previous equation produces the following Kelvin circulation theorem 



d 
It 



(u + R)-dx = 9(b C u {[z x (Vg)/77,u] • dx) 

c(u) J c(u) 



Modified PV advection. One way to obtain the modified advection equation for PV, is to take 
the exterior differential of the first equation in (3.12) (with a replaced by a in the second term) 
and obtain 




= 9d([z x (Vq)/n,u\ ■ dx). 



From this, and from the second equation in (3.12) one obtains 

d t q + u ■ Vg = 9-z ■ curl([z x (Vq)/rj, u]), 
by making use of the formulas 

— d ( — ] = q ■ dS and d(m ■ dx.) = (z ■ curl m) ■ dS 
for one-forms a and vector fields m on the two dimensional domain T>. 



Lagrange-d'Alembert variational principle. 

one obtains the RSW Lagrangian 



Making use of the Legendre transform (3.5), 



(u,r?) 



v 



+ u R 



D) 



dxdy . 
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From (3.6) we get the Lagrange-d'Alembert variational principle 

5 / £(u,n)dt + 9 / / rj([u, z X (Vq)/r)] ■ [u, v])dxdy, 
Jo Jo Jv 

for variations 5u = dtv + [u, v], Sr] = — div(^v), where v is and arbitrary time dependent vector 
field vanishing at t = 0, T, associated to a virtual displacement 5ip of the fluid configuration ip. 



3.2 3D rotating Boussinesq flows 

On the semidirect product of SDiff(IR 3 ) with functions b G J-"(IR 3 ), where SDiff(R 3 ) is the group of 
volume preserving diffeomorphisms of M 3 , the Lie-Poisson equations (3.1) become 

(Sh \ Sh Sh 
q x curl— ) - V— x V6 = 0, d t b + curl — ■ V6 = 0. (3.13) 
oq/ ob oq 

For the rotating 3D Boussinesq model, the Hamiltonian is 



h(q,b) = J Q| curP 1 q - R| 2 + bz*J d 3 x, 



where b is the scalar buoyancy and curl R = 2Q is the Coriolis parameter. The variational deriva- 
tives are ^ ^ 

curl — = curP 1 q — R = u, — = z , 
dq do 

and the vorticity is uj = curl u = q — curl R = q — 2f2. In this case, the Lie-Poisson equations 
(3.13) reduce to the Rotating 3D Boussinesq equations 

d t oj + curl ((w + 2O)xu)-zxVt = 0, d t b + u ■ Vb = , 

or, in velocity form 

<9fU + (curl u + 2fl) x u + zb = — Vp, 
where p is the pressure, found by enforcing incompressibility, div u = 0. 
The Lie-Poisson equations (3.13) admit the Casimir functions 

C$(q, b) = $(g, b) d 3 x, where q := V6 ■ q is the potential vorticity. 
Jv 

We shall now introduce the dissipative Casimir effect in the 3D Rotating Boussinesq system by 
using the Lie algebraic setting developed above. Choosing the enstrophy Casimir 

C(q,b) = ±J^q 2 d 3 x } (3.14) 

we have ^ ^ 

curl — - = curl(gV6), -— - = — div(gq) . 
dq ob 

Consequently, choosing the inner product 7 associated to a given positive self-adjoint differential 
operator A, the modified momentum (2.3) is 

q = q — 6 curl A [u, curl(gVfe)] = q — #L(curl(gV&) x u), 
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where the differential operator L verifies A = curl 1 L curl . The dissipative Casimir Lie-Poisson 
equations d t q + curl x curl Jj^ — V|| x V6 = 0, yields, in velocity form, the system 

d t u + curl (u + R - 0A[u, curl(gVfe)]) xu + bz = -Vp, d t b + u-Vb = 0. 

The first equation can be rewritten, in terms of L as 

d t u + (curl (u + R) - 6>L(curl(gV6) x u)) x u + fez = -Vp. (3.15) 

This system preserves the energy /i(q, b) of the 3D rotating Boussinesq fluid and dissipates the 
enstrophy as 



v 



9 || [curl(gVfe), u] || = — 9 ||curl(gVfe) x u 



i2 



(3.16) 



where we have used the identity (1.24). 



Remark 3.2. In the Casimir dissipation approach for the 3D Boussinesq equations, the vector 
curl(gVfe) in equations (3.15) and (3.16) now plays the role that the vorticity u> played for ideal 
incompressible 3D in equations (1.16) and (1.22). 



Kelvin circulation theorem. A convenient way to derive the Kelvin circulation theorem (3.3) 
is to rewrite the modification of equations (3.13) using duality pairing with one-forms. In this 
case, we have 

d t a + Cshdt — — db = —dp, 
6a ob 

in which a denotes the 1-form a = a ■ c?x or, more explicitly, 



d t a + CshO! = 9Ct 



Consequently, a direct computation yields 



Sh SC_ 

5a^ 5a 



+ —db - dp . 
db 



d 
It 



a = 9 



c(u) 



c(u) 



jC Sh 

da 



Sh sc_ 

Sol' 5a 



bd 



c{u) 



5h 
~8b 



consistently with (3.3), where we again recall that the Lie algebraic bracket is minus the Lie bracket 
of vector fields. For rotating 3D Boussinesq, upon using the enstrophy as a Casimir function and 
employing the differential operator A, one gets 



d 
dt 



j> (u + R) ■ rfx = 9 I £ u (A [u, curl(gVfe)] ■ dx) - j> bdz 

J c(u) J c(u) J c(u) 

= —9 (p u x L (curl(gVfe) x u) • dx — ® bdz, 

J c(u) J c(u) 

where we have used the expression for the Hamiltonian in terms of a and b, that is, 



h(a, b) 



v 



1, ^.9 , \ ^ $h 
- a - R + bz ft, — 
2 J da 



u, with ct = u + R . 



Setting 9 = recovers the usual Kelvin circulation theorem for rotating Boussinesq flows. 
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Lagrange-d'Alembert variational principle for Boussinesq flows. Proceeding exactly as 
in the preceding example for RSW and using the abstract formulation (3.6), one obtains the 
Lagrange-d'Alembert variational principle for Boussinesq flows, 



5 I £(u, b)dt + 9 [ I (curl(gV6) x u) • L(u x v)dxdy = 0, 
Jo Jo Jv 

for variations 5u = c^v + [u, v], 5b = — v ■ Vfo. Here v is an arbitrary time dependent vector field 
vanishing at t — 0, T, that is associated to a virtual displacement 5(p of the fluid configuration 
ip G SDiff(IR 3 ). Also, I is the Lagrangian of the 3D rotating Boussinesq fluid, given by 

£(u,b)= J Q|u| 2 + u-R-fej d 3 x. 



Remarkably, in the Casimir dissipation approach for 3D Boussinesq flows, the vector curl(gVfe) 
now plays the same role as the vorticity u> had played in the case of 3D ideal fluids. In particular, 
the vector field uj x u that appeared in the dissipative term for ideal fluids in the 3D Boussinesq 
case is replaced by the vector curl(gVfr) x u. 



4 Conclusion 

This paper has introduced a theory of selective decay by Casimir dissipation that applies widely 
in fluids and suggests a connection between selective decay and the design of numerical methods, 
since it recovers some of the previous methods, such as the anticipated vorticity method. 

Interestingly, while the entropy C 2 in (1.5) is the appropriate Casimir for selective decay in 2D 
incompressible flow, the appropriate Casimir for that role in 3D incompressible flow is the square 
of the helicity C in (1.13). 

Perhaps one reason for these associations with different Casimirs in 2D and 3D can be suggested 
by comparing the Lie-Poisson brackets in terms of the vorticity in 2D and 3D. In particular, these 
brackets may be written in terms of the Casimirs, the enstrophy C 2 in 2D and the helicity C in 3D. 
In fact, the Lie-Poisson brackets in 2D and 3D may both be written in terms of the corresponding 
Casimirs in the same triple product form Nambu [1973]. Namely, 



{/, h} 2 (u) 



r sc 2 


~5f 8h 


Iv 5u 


5uj' 5u 



dxdy and {/, h} 3 (oj) = [ curl — • curl -J- x curl ^— d 3 x . 

Jx> ouj ou: OUJ 



These formulas suggest the fundamental roles of the two different Casimirs, the enstrophy C 2 in 
2D for which 5C 2 /Su = u and the helicity C in 3D for which 5C/5ui = u. One could imagine using 
other constants of motion instead of the Casimirs, but these would result by Noether's theorem 
from particular symmetries and their use would be restricted to those subcases. 

The loss of both Kelvin's theorem and the Casimir conservation laws when modifying the 
equations to impose selective decay is consistent with the effects of viscosity, but of course the 
preservation of energy is not. This, however, is the essence of selective decay, which is seen in 
turbulence on times scales less than the time scale for viscous decay of energy. 

The phenomenon of selective decay is thought to arise from a multiscale interaction between 
disparate large and small scales in which the large scales are regarded as a type of coarse graining, or 
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perhaps the envelope, of the small scale motions, Mininni, Pouquet and Sullivan [2008]. This type 
of model is consistent with other recently developed multiscale turbulence models, particularly that 
in Holm and Tronci [2012], in which the large scale motion is regarded as a Lagrange coordinate 
for the small scale motion. 

Thus, a phenomenological point of view exists in which the effect of selective decay of Casimirs 
that depend on gradients of the velocity may be regarded usefully as providing a type of coarse- 
graining that is similar to the effect of a turbulent viscosity, but still conserves energy. In this 
regard, see also Vallis and Hua [1988]; Levy, Dubrulle and Chavanis [2006]. 

This interpretation of selective decay arising from coarse graining may also suggest a connection 
between selective decay and the design of numerical methods. In particular, one may wish to 
design numerical methods whose solutions imitate selective decay without adding too much linear 
viscosity. And from our discussion here it seems that the method of anticipated vorticity in 
Sadourny and Basdevant [1981, 1985] may be one of those numerical methods. 

The theoretical approach presented here may lead to classes of other numerical methods based 
on selective Casimir decay that would generalize and extend the applicability of the method of 
anticipated vorticity to other areas of fluid dynamics. A recent discussion of numerical methods 
based on anticipated vorticity for 2D flows is given in Chen, Gunzburger, and Ringler [2011a]. A 
"scale-aware" version of the anticipated vorticity method in 3D should also be readily obtainable, 
by allowing the value of its time-scale parameter 9 to depend on local mean properties of the 
solution. However, this is beyond the scope of the present work. 

Future applications of our approach here would naturally follow recent work in ocean dynam- 
ics that interprets selective decay as a mechanism for parameterizing the interactions between 
disparate scales, for example, between large coherent oceanic flows and the much smaller eddies 
directly, as done in Marshall and Adcroft [2010], instead of relying on the slower, indirect effects 
of viscosity. The Lie-Poisson structures in the selective decay models presented here suggest a 
framework in which numerical schemes for their computational simulations may be designed. This 
framework may turn out to be similar to the ideas in Salmon [2005] about using the triple product 
Nambu form to design simulation algorithms for certain grid structures used in oceanography. 
In addition, Voronoi mesh methods that were developed for anticipated vorticity dynamics may 
turn out to be naturally adaptable to the the selective decay models presented here. See. e.g., 
Chen, Gunzburger, and Ringler [2011a,b]. However, applications of these ideas for numerical sim- 
ulation algorithms within the Casimir dissipation modelling framework are beyond our present 
scope and will be pursued elsewhere. 

If selective Casimir decay or the anticipated vorticity method does mimic the effect of linear 
Navier-Stokes viscosity, but at a convenient time scale for numerical simulations, then one must be 
cautious about adding in other types of viscosity, such as the artificial viscosity used in Large Eddy 
Simulations to mimic the effects of turbulence. The introduction of other types of turbulent viscosi- 
ties such as Smagorinsky tensor diffusivity in addition to anticipated vorticity for selective decay 
would be seen as double counting. In contrast, the introduction of nondissipative regularization 
through nonlinear dispersion as a modification of the kinetic energy Hamiltonian, as accomplished 
in the Navier-Stokes alpha (NS-alpha) turbulence model Foias, Holm and Titi [2001], is comple- 
mentary to selective decay and cannot be regarded as double counting of dissipation. Moreover, 
the original derivation of the NS-alpha turbulence model by the method of Lagrangian averag- 
ing makes it clear that the model is already "scale-aware" because the length scale alpha in that 
derivation is the correlation length for Lagrangian fluctuations of the fluid parcel paths. 
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Analytical issues such as whether the selective decay equations admit existence of unique glob- 
ally strong solutions are not addressed here. However, in combining the NS-alpha model with the 
present selective decay approach, one would expect existence and uniqueness to persist, because 
this property is already possessed by the NS-alpha turbulence model. 

The Casimir dissipation approach for other fluid models such as magnetohydrodynamics (MHD), 
compressible fluids, two-layer quasi-geostrophic models parameterizing baroclinic instability, Maxwell- 
fluid and Maxwell- Vlasov plasmas, spin chains, solitons, complex fluids, turbulence models, etc., 
will be investigated elsewhere. 
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